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Abstract 

We study pair production of scalar top quarks (stop, ti) in e"'"e~ collisions with the 
subsequent decay of the top squarks into 6-quarks and charginos ii bxf. We simulate 
this process by using PYTHIA6.4 for the beam energy 2Ei, = ^/s = 350, 400, 500, 800, 1000 
GeV. A set of criteria for physical variables is proposed, which provides good separation 
of stop signal events from top quark pair production being the main background. These 
criteria allow us to reconstruct the mass of the top squark with an integrared luminosity 
of 1000 provided that the neutralino mass is known. 

1 Introduction. 

The scalar top quark, the bosonic partner of the top quark, has attracted much attention 
as it is expected to be the lightest colored supersymmetric (SUSY) particle, ti and t^R, 
the supersymmetric partners of the left-handed and right-handed top quarks, mix and the 
resulting two mass eigenstates ti and ^2, can have a large mass splitting. It is even possible 
that the lighter eigenstate ii could be lighter than the top quark itself [1]. 

Searches for top squarks which were performed at LEP and Tevatron [2], have shown that 
the mass of the stop is higher than 141 GeV for a mass difference between the stop and the 
lightest neutralino of about 50-70 GeV. These searches will continue at LHC and ILC [3], [4]. 

In the following we study the reaction ^ 

e+ + e" ii + ii . (1) 



^More details about this study can be found in [5]. The analogous analysis in the photon-photon channel 
was made in [6]. 



Among the possible ti-decay channels within the MSSM (see [7] for details), we focus on 
the decay ii bxf followed by the two-body chargino decay xf — ^ Xi^^; where one of the 
W's decays hadronically, W qiqj, and the other one leptonically, W jiu^ [8] ^. The final 
state of this signal process, shown in the left-hand plot of Fig.l, contains two 6-jets and two 
(or more) jets (originating from the decay of one W boson), a hard muon plus a neutrino 
(from the decay of the other W) and two neutralinos: 

e+e- t/i ^ hbxtxi ^ hhW^W-xlx\ ^ bbqiqjfiu.XiXi ■ (2) 

The main background process is top quark pair production with the subsequent decay t — >• 
bW^ (for W's we use the same decay channels as in the stop case): 

e+e- ~^ti^ bbW+W- bbqiq^fiu^ . (3) 

The only difference between the final states of stop and top production (shown in the right 
diagram of Fig.l) is that in stop pair production there are two neutralinos which are unde- 
tectable. Thus, both processes have the same signature: two 6-jets, two jets from W decay 
and a muon. In the present paper we consider only top pair production as background. The 




Figure 1: Left is the stop signal event diagram, Right is the top background diagram. 

analysis of the processes (2) and (3) is done on the basis of Monte Carlo samples of the cor- 
responding events generated by the two programs PYTHIA6.4 [10] and CIRCEl [11]. The 
program CIRCEl is used for the parameterization of the beam spectra involved in processes 
(1) and (3) to account for the effects of beamstrahlung. The energy of the beams was chosen 
as 2Eb = y/s = 350, 400, 500, 800, 1000 GeV. 

In Section 2 we give the set of MSSM parameters used in our study. 

In Section 3 we discuss some general characteristics of the signal process e+e^ tiii 
and the main background e~^e~ —>■ tt. First we show the result of taking into account beam 
interaction effects like beamstrahlung and others which are simulated with CIRCEl [11]. The 
distributions for a set of physical variables, which are measurable in the final state of the 
process (2), are given for the case of a beam energy y/s = 500 GeV and integrated luminosity 
of 1000 fb~^. These variables include hadronic jets from W-decay as well as 6-jets produced 
in stop decays. We always compare them in detail with those of top pair production. 

In Section 4 we propose three cuts which allow a good separation of signal stop events 
and background top events. The number of stop events that have passed these cuts and the 
values of the corresponding efficiencies are given. We show that the invariant mass of one 

^The process e+e^ ti +ti with the subsequent decay ti c^i was considered in [9]. 



6-jet and the other two non — 6-jets (from W decay) allows one to reconstruct the mass of the 
scalar top quark provided that the neutralino mass is known. The impact of the proposed 
cuts on the values of the cross sections of stop and top pair productions are shown together 
with the values of signal-to-background ratios (S/B). We have also considered a cut which 
allows to separate a signal muon from W decay and fake muons appearing due to decays of 
other particles in the same stop production event. 
Section 5 contains some conclusions. 



2 MSSM parameters and cross section. 

The scalar top quark system is described by the mass matrix (in the ii — in basis) [1], 



[12] 



M? M? 

tLL ii 

I? M? 

tRL tijfl 



with 

r2 71 /r2 



= ^0 + - \sin'Qw)cos2ml + Ml (5) 
Ml^^ = Ml + \sin^Qwcos2(3Ml + (6) 



The mass eigenvalues are given by 

1 



tl,2 

with the mixing angle 



cos9f = (9) 
'\M? P + (M^ -M? )2' ^ ^ 



Mf -M? 

sinOi = (10) 

In the following we will consider a particular choice of the MSSM parameters that are defined, 
in the notations of PYTHIA6.4, in the following way: 

M^ = 270 GeV; = 270 GeV; ^ = -500 GeV (top trilinear coupling); 
ton/3 = 5; = -370 GeV; Mi = 80 GeV; Ms = 160 GeV. 

Note that in PYTHIA6.4 corresponds to M^^ (left squark mass for the third gener- 
ation) and M^ corresponds to Mj^. These parameters give M^^ = 167.9 GeV, M^+ = 159.2 
GeV and M^o = 80.9 GeV. This value of Mj^ is rather close to the mass of the top quark 
Mtop = 170.9 ± 1.8 GeV [13]. Therefore one expects a rather large contribution from the top 
background, which means that the choice of this value of the stop mass makes the analysis 
most difficult. Finding a suitable set of cuts separating stop and top events is crucial. 

In general, the cross section for stop pair production at a fixed energy depends on the 
mass of the stop quark and the mixing angle 9i. Since the couplings of the to the left 
and right components of the stop are different, the cross sections depend significantly on the 



beam polarizations (see [7], [8], [14]). By choosing appropriate longitudinal beam polariza- 
tions it is possible to enhance the cross sections. For example, for an electron beam with 
90% left polarization the cross section would be larger than the unpolarized cross section by 
approximately 40%, for cos 9i = - 0.81 corresponding to the parameters given above. If in 
addition the positron beam has 60% right polarization, then the cross section is enhanced 
by approximately a factor of 2 compared to the unpolarized cross section. We note that a 
rather precise determination of the stop mixing angle 6j: is possible by measuring the left-right 
asymmetry. The cross section for top pair production has also a characteristic dependence on 
the beam polarizations [14]. For example, the polarization of both beams leads to an increase 
of the cross section by about a factor of 1.5. 



3 Distributions of kinematic variables in stop and top 
production. 

In this Section we present some plots of distributions for different physical variables based 
on 5 • 10^ stop pair production events generated by PYTHIA6.4 and CIRCEl weighted with 
the electron-positron luminosity. Analogous plots are also given for 10^ generated background 
top events. 

The ILC is a 200-500 GeV center-of-mass high luminosity (a peak luminosity of ~ 2 ■ 
10^'^cm~^s~^) linear electron-positron collider with a possible upgrade to 1 TeV in the second 
phase. According to [3] a total luminosity of 500 fb~^ is foreseen within the first four years 
of operation and 1000 fb~^ during the first phase of operation at 500 GeV. A first run at 
^/s = 500 GeV will get a first measurement of the particle masses to optimize the threshold 
scan [4]. 

Fig. 2 a) demonstrates the total energy spectrum of the electron and positron beams, 
which is expected at y/s = 500 GeV after taking into account beamstrahlung and other 
beam interaction effects (see, for instance [15]). Fig.2 b) shows the correlations of the beam 
fractions i/i = El^^^ (i=e+,e~) of the coUiding electron and positron beams. 



-500 



400 



300 



200 



100 





Entries 


1000000 










Mean 




494.3 










RMS 




10.2 










Underflow 















Overflow 















Integral 




1e+06 




































J 


...jj,... 


.... 


.... 




) 100 200 300 400 500 600 700 800 9001000 



Entries 


1000000 


Mean x 


0.9758 


Mean y 


0.9758 


RMSx 


0.03829 


RMSy 


0.03807 


Integral 


2.212e+05 





2498181280342 





221183248657 










0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



El, 



b) 



Vi 



Figure 2: a) total beam energy spectrum, b) beam fractions correlations 



To find the jets we used the subroutine PYCLUS of PYTHIA with the distance measure 
used in the "Durham algorithm". The parameters of this jetfinder are chosen such that the 
number of jets is exactly four. Technically, 6-jets are defined as jets that contain at least one 
B-hadron. Their decay may be identified by the presence of a secondary vertex [16]. 

Let us also mention that separate stop and top event samples were generated without and 
with cuts. This shows the effect of the cuts. 



3.1 Global jet variables. 



In this subsection we present some plots that were obtained without cuts. In addition to 
the two 6-jets there are also, according to the decay chain (2), two jets due to the decay of one 
W boson into two quarks W ^ Qi + qj (see Fig.l). It is well known that the characteristics 
of individual jets depend on the choice of the jet finder and its parameters. This dependence 
may result in some redistribution of energy between jets. To diminish these effects we use 
three physical variables which characterise the physical system composed of 6-jets and the 
jets produced by quarks in the W decay. The simulation has shown that the effect of energy 
redistribution between jets is partially compensated in multijet systems. 

For this reason we consider first the invariant mass of all four jets produced in the events 



which is the modulus of the vectorial sum of the 4-momenta Pj^^ ^ 



MinviAlljets) 



^j=l,2,3,4-P/et 



The distribution of this invariant mass is shown in Fig. 3. Plot a) shows the results for 
stop pair production while the plot b) is for top pair production. In Fig. 3 and the following 
figures the vertical axis shows the number of stop and top events that may be expected for the 
integrated luminosity of 1000 fb~^. Taking the integral of the distributions one can get the 
total number of events expected for the taken integral luminosity. These numbers are shown 
as "Integral" values in the Figures. One sees that for the chosen luminosity the number of 
the produced background top events (35 930) is about 15 times higher than the number of 
stop events (2373). It gives a signal-to-background ratio S/B = 0.066. 
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Figure 3: Distribution of events versus the reconstructed invariant mass of all jets Minv{AUjets). 
a) stop pair production, b) top pair production. 

Plot a) of Fig. 3 shows that the mean value of Minv{AUjets) is about three times lower 
in stop production than in top case, plot b) of Fig.3. Therefore, a rather soft cut hke 
Minv{AUjets) < 160 Gev, will strongly suppress the top background and save the bulk of 
signal stop events. Another useful variable that can also be used for the separation of the 
signal and the background is the "missing" mass (we use = 500 GeV) 



Mrr 



-j=l,2,3,4-t^j 



4i?j,, + £;^))2-(S,=i,2,3,4PL + P. 



(12) 



This variable takes into account the contribution of those particles that cannot be registered 
in the detector (neutrinos and neutralinos) . The distributions of this invariant "missing" mass 



3p = (Po, P), Po ^ Ep. We denote P^^^ as the S-momentum of the jet j: (P^), = (P-?);, {i = 1, 2, 3). The 
sum in (11) is over all four jets {j ~ 1, 2, 3, 4). 



are given in Fig.4. Plot a) shows the results for stop pair production, while the plot b) is for 
top pair production. As seen from these plots, the cut Mmiss > 250 GeV will also allow us to 
get rid of most of the background events. 
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Figure 4: Distribution of number of events versus the missing mass variable, a) stop pair production, 
b) top pair production. 



The third variable of interest is the invariant mass Minv{h-]ei.,J ET Sw) 
Minv{b — jet, JETSw) = Minv[b — jet + {all — non — b — jets)] = 

= \l {Pb-jet + Pall-non-b-jets)^ ; (13) 

which is constructed as the modulus of the vectorial sum of the 4-momentum Pb-jet of the 
6-jet, plus the total 4-momentum of "all-non-6-jets" system stemming from the W decay 
{Paii-non-b-jets = Pjetiw + -Pjet2M., as there are only two jets allowed to be produced in W 
decay). More precisely, if the signal event contains a /i~ as the signal muon, we have to take the 
6-jet (6-jet would be in the case of /i"*" signal muon) (see Fig.l). This is only possible if one can 
discriminate between b- and 6-jets experimentally. Methods of experimental determination of 
the charge of the 6-jet(6-jet) were developed in [18]. In this paper we do not use any b-tagging 
procedure. The PYTHIA information on the quark flavour is taken for choosing the b and b 
jets. 

In the top case the invariant mass Minv{b, 2 quarksw) of the system composed of a 6-quark 
and two quarks from W decay should reproduce the mass of their parent top quark (see Fig.l). 
The distributions of events dN^""^"^^ / dMinv/^ GeV expected in each bin of 5 GeV versus the 
invariant mass Minv{b,2 quarksw) of the parent three quarks as well as the invariant mass 
of jets produced by these quarks, i.e. Minvip — jet, JETSw), are shown for quark and jets 
levels in the plots a) and b) of Fig.5, respectively, for integrated luminosity of 1000 fb~^. 

The distributions of Fig.5 show that the peak positions at quark level (plot a) ) as well as 
at jet level (plot b) ), practically coincide to a good accuracy as well as with the input value 
of the top quark mass Mtop = 170.9(±1.8) GeV. 

It is also seen from the plot b) of Fig.5 that the quark hadronisation into jets leads to a 
broadening of very small tails which are seen in the invariant mass distribution at quark level 
(plot a) of Fig.5). The right tail is a bit lower and longer than the left one. One may say that 
the peak picture at jet level still looks more or less symmetric. The main message from these 
plots is that the appearance of tails due to quark fragmentation into jets does not change the 
position of the distribution peak, which allows us to reconstruct the input top mass both at 
quark and jet level. 
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Figure 5: The top pair production spectra of the invariant masses Minv(b-jet,JETSw) of the "b- 
jet+(aU- non-b-jets)" system, a) at quark level, b) at jet level. 



4 Cuts for the reconstruction of the stop mass. 

Based on the results of the previous Section we will use the following three cuts to 
separate the signal and background: 

• there must be at least two 6-jets in an event: 

Nb-jets > 2; (14) 

• the invariant missing mass must be larger than 250 GeV: 

Mrniss > 250 GeV- (15) 

• the invariant mass of all jets must be smaller than 160 GeV: 

Minv{AUjets) < 160 GeV. (16) 

They reduce the number of top background events from 3.5*10^ to 12 and leave 1806 stop 
events. So, the cuts improve the signal-to-background ratio from S/B = 0.066 to S/B ^ 143 
losing about 24% of the signal stop events. The efficiency values for the cuts (14)-(16) are 
calculated for ^/s = 500 GeV. We define them as the summary efficiencies. It means that if 
El is the efficiency of the first cut (14), eu is the efficiency of applying the first cut (14) and 
then second cut (15). Analogously, £123 is the efficiency of the successive application of the 
cuts (14), (15) and (16). The following results are obtained: 



For SIGNAL STOP events: ei = 0.84; = 0.78; £123 = 0.76; 

For BACKGROUND TOP events : Ei = 0.94; £12 = 0.001; £123 = 3.5 ■ 10"^ 

The distribution of the invariant mass of the "6-jet+(all-non-6-jets)" system in the case 
of stop pair production is shown in Fig.6. The plots of this figure are analogous to the 
corresponding plots of Fig. 5 using only those stop events that have passed the cuts (14)-(16). 



In the stop case, one has to take into account that one of the stops decays into three jets 
plus a neutralino Xv Therefore, the right edge of the peak of invariant mass distribution of 
the " b-}et+{J ET Sw*)" system corresponds to the mass difference Mj^ — M^o. The position 
of the right edge of the peak Mj„„(6-jet, JETSw*) ~ 85 GeV (plot b) of Fig.6) coincides 
with the peak position seen at quark level (plot a) of Fig.6). Let us note that the observed 
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Figure 6: The spectra of the stop signal events after cuts versus the invariant mass Minv (b-jet, 
JETSw* )■' a) o-t quark level, b) at jet level. W* means virtual nature of W boson. 



stability of the peak position in both plots of Fig. 6 is due to the rather moderate loss of the 
number of events that form this peak (they change from ^ 550 to ~ 280) while passing from 
quark level to jet level. Adding the mass of the neutralino M^o = 80.9 GeV one gets the 
reconstructed stop mass M~^'^° k, 166 GeV which reproduces well (within bin width of 5 GeV) 
the input value Mj^ = 167.9 GeV. 

The simulation has shown that the 12 background events are mostly distributed in the 
region 40 < Mi„^(6-jet, JETSw)< 140 GeV, which is 20 times wider than the 5 GeV width of 
the peak which contains about 280 signal stop events left after the cuts. Therefore, we expect 
that in future measurements the contribution of a few remaining top background events will 
not influence the position of this peak which allows one to reconstruct the input value of the 
stop mass by adding the mass of the neutralino. 

The origin of the left tail of the distribution shown in the plot a) of Fig. 6 can be clarified 
by the results of the stop mass reconstruction by calculating its invariant mass at quark level 
Minv{b, 2 quarksw*, Xi) the modulus of the sum of the 4-momenta of all three quarks and 
the neutralino (see Fig.l) produced in stop decay. These results are given in the plot a) of 
Fig. 7 which shows a very precise reconstruction of the input stop mass at quark level withing 
the 5 GeV width of a peak. 
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Figure 7: The spectra of the stop signal events after cuts versus the invariant mass Mi, 
JETSw* ,Xi) ^) o-t quark level b) at jet level. 



, (b-jet. 



Comparing the plot a) of Fig. 6 with the plot a) of Fig. 7 one can conclude that the long left 
tail as well as the very small right tail in the distribution of Mj„t,(6, 2 quarksw*) are caused 
by the fact that the neutralino 4-momentum is not included in the "6-jet+(all-non-6-jets)" 



system. Analogously, comparing plots b) and a) of Fig. 6 one may estimate the contribution 
of the hadronisation effect (in addition to the contribution of the neutralino momentum) to 
the left tail of the distribution b) in Fig.6. 

The influence of the effect of the hadronisation of the 6-quarks and of the quarks from 
W decay into jets is shown in the plot b) of Fig. 7. It is seen that the hadronisation of 
quarks into jets does not change the position of the stop mass peak, which still includes the 
input value M^^ = 167.9 GeV, but more or less symmetrical and rather suppressed short tails 
appear around the peak position. The shape of the peak in the stop plot b) Fig. 7 looks 
very similar to the shape of the peak in the top plot b) of Fig.5, which demonstrates the 
stability of the reconstructed top mass peak position after taking into account the effect of 
quark fragmentation. 



4.1 Separation of signal and fake muons. 

To select the signal stop pair production events shown in the left plot of Fig.l, one has 
to identify the muon from the W decay. The distribution of the energy of the signal muons 
Esig-mu, obtained from the sample of events generated without any cuts, is shown in Fig. 8 a). 
There are, however, also muons in the event coming from leptonic and semileptonic decays of 
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Figure 8: a) Energy distribution of signal muons. b) Energy distributions of fake muons. 



hadrons. Fig. 8 b) shows the energy Edec-mu of these muons stemming from hadron decays 
within the detector volume (for which we took the size from [3] , [4] ) . It can be seen that the 
decay muons have a rather small energy Edec-mu- Its mean value is about 0.81 GeV. The 
analogous spectrum for the signal muons in Fig. 8 a) shows that the signal muons have a much 
higher energy. The mean value Esig-mu is about 50 times higher than the mean value of the 
energy of the decay muons. One can cut off most low-energy decay muons rejecting those 
with Emu < 4 GeV. It leads to a loss of about 15 signal events as seen from the plot a) of 
Fig. 8 (the bin in this plot is 2 GeV). 

We have also studied another way to select the signal muon from W decay. If the axes of 
all four jets in the event are known, then in general the signal muon has the largest transverse 
momentum with respect to any of these jet axes. 



4.2 Signal and background cross sections. 

We give in Table 1 (without cuts) and Table 2 (after the cuts) the cross sections and the 
numbers of events for stop pair production and top pair production for five energies i/i. To 



estimate the rates of stop production at different energies we take the universal value of the 
luminosity of 1000 fb~^ for all energies. 

It is seen that the fixed cuts (14)-(16) lead to a strong (more than 3 orders of magnitude) 
suppression of background top contribution and a moderate loss of signal stop events produced 
in the energy range 400 < ^/s < 800 GeV. 

For the stop mass chosen, the largest number of signal events is expected at ^/s = 500 
GeV. Let us note that, according to [18], a 50% efficiency of the separation of b jets and 
80% of the corresponding purity can be expected. It means that to get 1806 reconstructed 
signal stop events wc will need about 2.5 times higher statistics than that provided by the 
luminosity 1000 fb"^ at = 500 GeV. 

It is worth noting that, as seen from Table 2, even with the use of the fixed parameters of 
the cuts (14)-(16) the number of signal stop events that can pass these cuts grows rapidly with 
the energy in the region 400 < y/s < 500 GeV where mass measurements will be done in the 
first phase. These measurements may allow to enlarge the sample of collected signal events 
and to perform a precise measurement of the stop mass. The region 500 < y/s < 800 GeV, 
which be available in the second phase of ILG operation, allows to gain a much higher stop 
statistics, as seen in Table 2. A complete analysis based on adjusting the parameters of the 
selection cuts (15)-(16) for each of the above energy intervals will be presented in forthcoming 
papers. 

Table 1: The cross sections and the number of events for stop and top pair production before 
cuts. 
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Table 2: The same as in Table 1 but after cuts. 
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5 Conclusion. 

Wc have studied stop pair production in electron-positron collisions for a stop mass of 
167.9 GeV within the framework of the MSSM for the total energies 2Eb = \fs = 350, 400, 
500, 800, 1000 GeV. We assume that the stop quark decays dominantly into a chargino and a 
b quark, ii bxt, and the chargino decays into a neutralino and a W boson, Xi ~^ Xi^^? 
where the W boson is virtual. One of the two W's decays hadronically, W'^ — > qq, the other 
one decays leptonically, W' — > 



We have performed a detailed study based on a Monte Carlo simulation with the program 
PYTHIA6.4 for ^/s — 500 GeV (at this energy we expect the highest number of the signal 
events for the chosen stop mass Mf^ = 167.9 GeV and cuts) and the luminosity 1000 fb~^. 
The program ClRCEl is used to get the spectra of electron (positron) beams taking into 
account the effects of beamstrahlung. PYTH1A6.4 is used to simulate stop pair production 
and decay as well as top pair production being the main background. 

Three cuts (14)-(16) for the signal stop events have been proposed to separate signal stop 
events and top background events. For i/s = 500 GeV and the luminosity 1000 fb~^ they 
give 1806 signal stop events with 12 top background events. This is different from the more 
complicated situation in stop pair production at LHC (sec, for instance, [17]). 

We have shown that the determination of the peak position of the distribution of the 
invariant mass Mj„^(6-jet, JETSw) of "6-jet+(all-non-6-jets)" system allows one to measure 
the mass of the stop quark with a good accuracy based on the statistics corresponding to the 
luminosity 1000 fb~^. For this the mass of Xi has to be known. 

As seen from the Table 2 the measurements at other energies in the regions 400 < < 500 
GeV and 500 < ^/s < 800 GeV may allow to enlarge substancially the number of selected 
signal stop events and to perform a precise measurement of the mass of the scalar top quark. 

In conclusion we can say that the e^e~ channel is well suited for the study of stop pair 
production at ILC. 
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